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Abstract 

The effect of the nucleon recoil for antikaon-deuteron scattering is investigated in the framework 
of effective field theory. In particular, we concentrate on the calculation of the nucleon recoil effect 
for the double-scattering process. It is shown that the leading correction to the static term that 
emerges at order ^^/^ with ^ = MkIttin vanishes due to a complete cancellation of individually 
large contributions. The resulting recoil effect in this process is found to be of order of 10-15% 
as compared to the static term. We also briefly discuss the application of the method in the 
calculations of the multiple-scattering diagrams. 
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1 Introduction 



Recent combined analyses of pionic hydrogen and pionic deuterium have revealed a significant progress 
towards a precise extraction of the pion-nucleon interaction parameters at threshold, see Ref. [Ij and 
also [2j ,3J for review articles. One might expect that the properties of the KN interaction can be 
investigated in a similar way using a combined analysis of kaonic hydrogen and kaonic deuterium. The 
data on this sort of bound systems are provided by the ongoing experiment of the DEAR/SIDDHARTA 
collaboration, which aims to measure the Is energy level shift and width of kaonic hydrogen and kaonic 
deuterium, eventually with an accuracy of several eV [3| . Contrary to the vrA^ case, the KN scattering 
length, however, is known to be large (of order of 1 fm) and strongly absorptive. This means that in 
this case one will need to measure 4 independent quantities (the real and imaginary parts of the S- 
wave KN scattering lengths bo, bi), whereas for vrA^ scattering one has to deal with two real scattering 
lengths. Consequently, the role of the deuterium measurements is rather different in these two cases. 
While the pion-nucleon scattering lengths can, at least in principle, be determined from the pionic 
hydrogen alone with the data from pionic deuterium serving as an additional consistency check, both 
data from the kaonic hydrogen and kaonic deuterium are needed to determine the values of ^i- 

The analysis of the experimental data proceeds as follows. From the measured values of the energy 
shift and width of the ground state in the kaonic hydrogen (kaonic deuterium) one first extracts the 
threshold scattering amplitudes of K'p {K~d) scattering by using the DGBT-type formulae [5] at 
next-to-leading order in isospin breaking (see, e.g., [21 El [7]) 

/S.Eis — i—^ = — 2a^/x^ap(l — 2^a(ln Q — l)ap) , 

/\Ei-i^ = -2aVd^;^d(l-2/^<i«(lna-l)A^rf), (1) 

where ^ {fid) stands for the reduced mass of the {Kd) system and Op {Aj^^) refers to the pertinent 
threshold amplitude^. Here and in what follows, we adopt the notation of Ref. [2]. 

In the next step, one has to express the quantities Op and Af^^ in terms of the S-wave KN scattering 
lengths bo, bi which are defined in the isospin-symmetric limit at q = and rrid — niu = 0. The isospin 
structure of the KN scattering amplitude in this world is proportional to bo + biTK-^N and the relation 
to the scattering lengths aQ,ai corresponding to the total isospin I = 0, 1 is given by = bo — 36i, 
oi = ^0 + bi- Our ultimate goal is to extract the precise values of bo,bi from the experiment which 
are then to be confronted with the theoretical predictions obtained in the unitarized ChPT [9]. This 
would provide a beautiful test of our knowledge about the SU{3) meson-baryon dynamics at low 
energy. Eventually, it will be very interesting to compare the experimental results with the scattering 
lengths directly extracted from lattice QCD with the use of the recently proposed method [lOj . 

In the isospin-symmetric world, Op = bo — bi. This relation is strongly modified when isospin is 
broken due to the unitary cusp effect [6l[Tl]. However, the corrections due to the cusp effect can be 
resummed to all orders in perturbation theory leading to the modified expression for Op which is still 
given in terms of bo, bi and the physical masses of the kaons and nucleons [U [11]. Thus, the presence 
of the cusp effect in ap does not affect the accuracy of the extraction of ^o^^i from the data. 

^In the case of kaonic atoms, higher-order Coulomb corrections may turn out to be not completely negligible numer- 
ically as shown, e.g., in Ref. [8] through the exact solution of the Schrodinger equation for the bound state. This issue 
is, however, relatively easy to cure since the large contribution comes from an iteration of a particular diagram to all 
orders. Replacing the factor 1 — 2fia{ln a — l)ap by (1 + 2^a{\n a — l)ap)~^ already captures the bulk of the effect. We 
shall not further elaborate on this issue. 
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The situation is rather different for the kaonic deuterium, even in the absence of isospin breaking. 
In order to perform the analysis of the data, the quantity Aj^^^ should be expressed in terms of 60, bi 
in some kind of multiple-scattering expansion. The procedure is plagued by our poor knowledge of 
the 3-body Kd interactions. It is not clear a priori whether the corresponding uncertainties are large 
enough to preclude one from being able to extract the values of bo, bi from experimental data if such 
precise data would be available. A systematic model-independent analysis of the uncertainties in the 
3-body sector is therefore needed in order to properly analyze the kaonic deuterium data, which will 
be provided by SIDDHARTA collaboration in the near future. 

The aim of the present work is to address the recoil effect, which is one of the major sources of the 
theoretical uncertainty in the study of Kd scattering at low energy. Our paper is organized as follows. 
In section[2]we briefly review the theoretical framework, which will be used to systematically investigate 
recoil effect in the Kd scattering. The background information about the recoil effect is given in 
section [3l Section H] is devoted to the discussion of the recoil effect in the double-scattering diagram and 
contains both the analytic calculations and numerical results. Recoil effect in the multiple-scattering 
diagrams is briefly discussed in section [5l Finally, our conclusions are presented in section [H 

2 The framework 

During the last few decades, the Kd scattering at low energy has been addressed on numerous occasions 
within the framework of Faddeev equations, see e.g. Refs.[T2l [T3\ [T^ [T5] for some of those works. It 
should, however, be pointed out that the results of these beautiful calculations are of no direct use 
in the analysis of the SIDDHARTA data since these calculations do not provide an explicit relation 
between Aj^^ and bo, bi, which is needed for the analysis. 

Contrary to the brute-force Faddeev calculations, the multiple-scattering series for the Kd scattering 
length, see e. g. [3 \W[ I17j . can be utilized for the analysis of the data. However, as it is well known, 
the KN scattering lengths are large due to the presence of the subthreshold A(1405) resonance, and 
the multiple-scattering series does not converge. The resummation of the series can be carried out 
analytically using the so-called fixed center approximation (FCA) in which nucleons are considered 
as static sources, i.e. mjsf 00. However, since Mx/rriM — 0.5, one may a priori expect large 
corrections to the static limit. To the best of our knowledge, no systematic method to evaluate these 
recoil corrections exists in the case of a large scattering length, where the multiple-scattering series 
should be resummed. The present paper aims at the formulation of the framework, which is capable 
to address this problem. 

A comparison of the exact numerical solutions of the Faddeev equations for potential models with 
the multiple-scattering series resummed in the FCA reveals a pretty good agreement in most cases, 
see |18] and references therein. This observation provides motivation for our approach and serves as a 
starting point since it indicates that the recoil corrections could be rather small even at M/^/mjy — 0.5 
and might be amenable to the perturbative treatment. Note that the perturbative treatment of the 
recoil corrections becomes indispensable if the resummed multiple-scattering series is considered. The 
potential model is a useful testing ground for our framework, since it allows to examine the convergence 
of the perturbative series towards the exactly known result. 

As pointed out in Ref. [7J, the non-relativistic effective field theory (EFT) provides an ideal tool to 
explore the multiple-scattering expansion. In particular, one readily reproduces the known result [16^ 
[T7] for static nucleons, and the inclusion of the recoil effect is straightforward (at least, formally). 
Prior to going to the details of the calculation we find it appropriate to discuss the essential features 
of the EFT approach to the problem of interest. 

i) The usefulness of the multiple-scattering expansion for the Kd scattering is due to the ap- 
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pearance of two distinct momentum scales. Whereas the A^A^ interactions and 3-body KNN 
interactions are mediated at large distances by 1-pion exchange, the dominant long-distance con- 
tribution to the scattering is governed by 2-pion exchange. For this reason, we may treat 
KN interactions as point-like, whereas A^A^ and KNN interactions are described by non-local 
potentials. 

ii) The fact that A(1405) resonance is located close to the KN threshold can, potentially, lead to 
enhancement of the non-local contributions in the KN scattering amplitudes. One may expect, 
however, that the non-local effects could be still taken into account perturbatively, by using the 
effective-range expansion of the KN amplitude. This assumption can be checked a posteriori 
by explicit calculations. 

iii) The expansion parameter associated with the ^A^ interaction is given by b ■ (l/r) ~ 1, where h 
denotes the S-wave KN scattering length and (l/r) ~ 0.5 fm~^ is the expectation value of the 
operator r~^ between the deuteron wave functions. Since the expansion parameter is large, the 
multiple-scattering series does not converge and should be resummed to all orders. 

iv) At the energies which are relevant for the problem in question, effects due to explicit cre- 
ation/annihilation of the particles can be neglected. These processes are taken into account 
implicitly through various low-energy constants and non-local interactions present in the La- 
grangian. The non-relativistic EFT, which we use, conserves the hadron number explicitly. 

v) In addition, we expect that there is no need to include hyperonic channels KNN — ttYN with 
y = A, S explicitly in the non-relativistic EFT framework. This conjecture is backed by the 
numerical coupled-channel Faddeev calculations [121 13 13 HI] • 

vi) The (non-local) A^A^ potential can be directly imported from chiral effective field theories, see, 
e.g., [19 ^ 120 ^ [2T] for recent review articles. The A^A^ and KN sectors of the non-relativistic EFT 
do not talk to each other. However, in this work for demonstrative purposes we shall use NN 
potential in the separable form. 

vii) Inclusion of the three-body force is required in EFT for the consistency reasons alone. In 
order to estimate the numerical strength of the three-body force notice that the total two- 
nucleon absorption rate in K~d scattering amounts for (1.22 it 0.09)% [22]. Assuming that the 
dispersive and absorptive parts of the 3-body force have the same order of magnitude, we may 
expect effects of the 3-body contribution to Aj^^ at a few percent level. This is definitely beyond 
the theoretical accuracy of the present calculation. For comparison, note that the imaginary part 
of the nd scattering length, which is dominated by the contribution from vrd NN breakup 
reaction, amounts more than 20% of the real part |23j . 

viii) As we shall see below, the use of the EFT framework enables one to systematically obtain the 
expansion of Aj^^ in powers of the inverse nucleon mass. The leading-order term corresponds 
to the static nucleon limit with — > oo. The key assumption is that this expansion is 
perturbative, even if the multiple-scattering series is not. 

ix) In order to simplify the bookkeeping, we neglect relativistic effects for the time being. Stated 
differently, we consider EFT of the underlying fundamental theory which is assumed to be 
the non-relativistic potential model. All general conclusions should remain in place for such a 
model. A systematic inclusion of the relativistic corrections will be discussed elsewhere. In the 
non-relativistic case at hand one may introduce the dimensionless parameter ^ = Mx/mN and 
consider the expansion of the amplitude ^/f^ in this parameter. Both integer and half-integer 
powers of ^ appear in this expansion, emerging from different momentum scales. 
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Figure 1: Diagrams corresponding to the meson double scattering on the nucleons in the deuteron 



x) For the time being, we completely neglect isospin-breaking effects in Kd scattering. They will 
be considered in a later publication. 

Equipped with this general clue, we now consider the perturbative calculation of the recoil correc- 
tions in detail. 

3 Some known facts about nucleon recoil 

Historically, the appearance of sizable cancellations in the case of the vrd scattering length was first 
pointed out in the papers by Kolybasov et al. and, independently, by Faldt [25] who argued that 
the naive static term provides a good approximation of rescattering effects. Recently, the role of the 
recoil effects was investigated for vrd scattering and for the reaction — > ir^nn within EFT |26l I27j . 
In particular, it was shown that the 3-body singularity that occurs in the ttNN intermediate state is 
the origin of non-analytic corrections to the static term of order y^M^Jyrn/v . It was also found that the 
importance of the recoil effect is directly connected with Pauli selection rules for the intermediate A'^A^ 
state. In the case when the S-wave NN interaction is forbidden by the Pauli principle, the leading 
MT^/mjsf correction cancels in the diagrams of double vrA^ scattering, diagrams a) and b) in Fig. [U 
thus resulting in the small net contribution. The situation is different in case when the S-wave A^A^ 
interaction is allowed, see diagram c) in Fig.[Tl First, the leading MT^/rriN correction does not vanish 
here (the contributions from diagrams a) and b) appear with the same sign). For this reason, it was 
concluded in [26^ [27| that recoil effects should be significant in this case. Secondly, the additional 
diagram with the S-wave NN interaction in the intermediate state must also be taken into account 
in the calculation. In the process ^d ir^nn, the intermediate A^A^ interaction appears to be in 
the ^Sq partial wave. Since A^A^ interaction in the ^Sq partial wave differs considerably to the one 
in the ^Si channel there is no a priori reason to expect any kind of cancellation between the recoil 
correction from diagrams similar to a) and b) and the contribution of diagram c). However, for nd 
and Kd scattering the intermediate and final A^A^ interaction occurs in the same channel {^Si —^Di). 
Therefore, a combined consideration of all diagrams of Fig.l is needed to draw a definite conclusion 
about the recoil effect!. Note that the Pauli-allowed recoil correction for vrd scattering gets multiplied 
with the isoscalar ttN scattering length squared which renders this effect negligible. This is, however, 
not the case for Kd scattering where both the isoscalar and isovector interactions are of a similar (and 
large) size. 

^We neglect the D-wave contribution in the foUowing. This does not affect the structure of the recoil contribution, 
which is of our primary interest here. 
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4 Kd scattering: Recoil effect in the double-scattering process 
4.1 The method 

Despite the well-known fact that the multiple-scattering series for the Kd scattering requires resum- 
mation to all orders, at the first stage of our investigation we concentrate on the double-scattering 
diagram. The aim is to demonstrate the technique used to obtain a systematic expansion of any 
Feynman diagram in (generally non-integer) powers of the parameter ^. To this end, we apply the 
perturbative uniform expansion method proposed in Ref. [28] . The same results can be also obtained 
using the threshold expansion method for the Feynman diagrams developed in Ref. p9]. The advantage 
of the scheme of Ref. ^28j is, however, that it is not tied to a particular (dimensional) regularization. 

Application of this method to carry out the integral in an arbitrary one-loop Feynman diagramU 
can be summarized as follows: 

i) Identify the relevant momentum scales and separate the range of integration into the regions 
according to these scales. Suppose, for instance, that there are two distinct scales: the small 
scale r] and the large scale A, 77/A <C 1. Then, the integrand f{r],q,A) with q referring to the 
integration momentum has three relevant regimes: 

1. The low- momentum regime with r/ ~ g ^ A, 

2. The high- momentum regime with 77 ^ (7 ~ A, 

3. The intermediate regime with rj <^ q <^ A. 

ii) In each region, perform the Taylor expansion of the integrand /(r/,g. A) in the corresponding 
small parameters. 

iii) Then, the original integrand fulfills the equality 

fiv, q, A) = Mr,, q, A) + Mr,, q. A) - Mr,, q, A) (2) 

at the given order in r]/A. Here, f\, /h, /i, refer to the Taylor-expanded integrand / in the low-, 
high- and intermediate-momentum regimes, respectively. In particular, the integrand /i, which 
represents the intermediate region, contains infrared- and ultraviolet-divergent terms necessary 
to make the integrals over the functions fi and fi finite and thus plays the role of a regulator. 

iv) Finally, integrate the functions /i,/h,/i over the whole momentum range. Since the above 
combination of these functions reproduces the original function f{r], q, A) at the given order in 
ri/A, the same combination of integrals will yield the correct result for the original integral up 
to the same order. 

Let us now apply this method to the double KN scattering diagrams of Fig. [TJ The corresponding 
contribution to the Kd scattering length reads, see Ref. ^26j for more details: 

Ap/- = ^^^^ {Ra + R, + Rc) , (3) 



''The method is, however, not restricted to one-loop diagrams. 
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where the quantities Ra, Rb and Rc are defined as follows: 



Ra 

Rb 
Rc. 



2Mk 
2Mk 



(27r)6 



P+2 



2Mk 



(27r)e 



rriN 
12 



4m7v 



-1 


ri2i 




2/i 



2M/^ 



+ 



p2 + 7^ 
rriN 



+ 



12 



4m7v 



4m2^M;^ 
12 



(27r)£ 



2Mi^ 



+ + 

rriN 



12 



4?TT,7V 



-1 



^(p+- ) ^'(q+- )Mw^(p,q,£;(l)) 



q2 + 72 12 
+ + 



(4) 



where ^ denotes the deuteron wave function, 72 
and the energy £^(1) is defined as 



2Mk m N 4m n 

mN£d with Ed being the deuteron binding energy 



12 



12 



2Mk 'irriN 



(5) 



Further, M^^ denotes the A^A^ amplitude. Its normalization is chosen such that in the CM frame, 
the amplitude Mjsfjsf evaluated on the energy shell is related to the S-wave scattering phase shift 5{k) 
through Mj^N{k,k,k'^ /rriN) = 16TTmi\f{k cot6{k) — ik)~^ . In the quantity Ra, renormalization of the 
one-loop KN scattering amplitude is carried out by performing the subtraction at threshold. At this 
order, this prescription yields the same result as dimensional regularization, but has the advantage 
that it is not tied to a particular regularization scheme. 

In the static limit with ^ = Mk/thn — > 0, only the amplitude Rb survives. In the vicinity of the 
static limit, each of the amplitudes Ri, i = a,b,c can be expanded in half-integer powers of ^ 



= i?f ^ ^1/2 ^ ^ p(2) ^ ^3/2 ^(3) ^ _ _ _ _ 



(6) 



Our aim is to perform this expansion systematically using the perturbative uniform expansion method 
[28j . In order to do this, it is convenient to rewrite the amplitudes Ri by showing the explicit depen- 
dence on the parameter ^ in the integrands 



Ra + Rb + Rc — &o 



where 



'o(i) 



P + 



Inn 



d^p dH 
(27r)6 



C fd^pd^qdH 



±1- ( P + -) ^ ( p 



l' + e(2(p' + 7')+lV2) 1^(1+0 



TUN 



I 



^(P+^) *(q+^ )M^^(p,q,^(l)) 



(27r)9 
(fpdH 

(2^)6 



[12 + e (2(p2 + 72) + 12/2)] [12 + ^ (2(q2 + f ) + l2/2)] ' 



*(p + ^)*(p-^)^, 



'st 



(7) 

, (8) 

(9) 

(10) 



Here, 1st corresponds to the FCA result (obtained in the static limit). The recoil corrections corre- 
sponding to the Pauli-allowed (forbidden) S-wave A^A^ intermediate state in the diagrams a) and b) of 
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Fig. [Dare given by the integrals Iq and A/st (/i and A/st) in Eqs. ([8]) and (fTO]) . For the Pauli- allowed 
A^A^ state there is also a contribution from the diagram c) given by the integral Inn- 

We define the recoil corrections as AIi = — /i + A/st and AIq = /q + /tvat + A/jt for isovector 
and isoscalar ^A'^ interactions, respectively. Below, we demonstrate how these corrections can be 
evaluated by using the method of Ref. [28] . 

Consider, for instance, the integral Ii, which can be rewritten in the following form 

( ^<iW./(P,.) ;(p.i) = ^(p.l) „ '7i-t:„„ . (U) 



' l + U (2vr)6 P ' •'^^'^ '12 + ^62 + ^12/2 

Here, = 2(p2 + 72) and the quantity <I> denotes the following combination of the wave functions 

$(p,l) = *2|^p + ^^_^j^p+l^^(^p_l^ (12) 

There are three relevant momentum regimes in the integral given by Eq. (jlip : 
1. The low-l regime. 

In this case, the involved momenta scale as follows: 

l2 p2 

r =^ 1~\/^P, P~(lA)wf, (13) 



2Mk 2m at 

where (. . .)^f denotes averaging over the deuteron wave functions. In this regime, ^/^b ~ 1 and 
1 <C 6 ~ p. Consequently, there are two different expansion parameters: ^ and I'^/b'^. Note, however, 
that the term ^b'^ occurring in the propagator of Eq. (jlip is of the same order as P and thus should 
be kept unexpanded. Expanding the wave functions in Eq. (jl2p in powers of the momentum 1 and 
using <I>(p, 0) = we get 

ci>(p,i) = ^ia,vivfc^(p,i) +.... (14) 
^ 1=0 

Note that we only keep terms even in 1 in the expansion of $(p, 1) since all odd terms will vanish after 
performing angular integration in Eq. (jlip . Expanding the integrand /(p, 1) and averaging over the 
directions according to Ijlj — > 5ijl'^/3 leads to 

/l(p, 1) = cl>2(p) (^-62 + 12/2 + J^f^) + ■ • • With <I.2(p) ^ ^V2<I>(p, 1) 1^^^ , (15) 

where the ellipses refer to terms of a higher order in ^. We attach the subscript "1" to the expanded 
integrand in order to signify the low-l regime. 

2. The high-l regime. 

In this case, the momenta 1 and p scale as follows: 

l~p-(l/r)wf =^ 7e^«l~6. (16) 
This implies that the function /(p, 1) can be safely expanded in powers of 

A,p,)..,,„(z^,,^),.... a„ 

We attach the subscript "h" to the expanded integrand in the high-l regime. 
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3. The intermediate regime. 

This case corresponds to the scahng ^/^b ^1^6. The expanded integrand /i in this regime can be 
obtained by using the heavy-baryon expansion to the function fi or the low-momentum expansion to 
the function /h- Both expansions lead, of course, to the same result 



/i(p, 1) = cl>2(p) ( -b' + 1^2 + C j2 ) + . . . . (18) 

Adding up the contributions from all three different regimes we obtain 

It is now easy to identify the powers of ^ which emerge after carrying out the integration. The first 
two terms are polynomials in ^. Recalling that the whole integral is multiplied by a factor C/(l + it 
is seen that these terms start to contribute at 0{S,) and 0{^^), respectively. Rescaling the integration 
momentum 1 — > 1 in the third term, one sees that this term contributes at order 

It is easy to verify that the method of Ref. |28j indeed leads to a systematic expansion in ^ by 
noting that the neglected terms A/ = / — (/i + /h — fi), 

yield contributions to Ii in Eq. ()lip of order of 0{^^) after evaluating the corresponding integrals. At 
first sight, one might expect that rescaling 1 — > -y/f 1 in the second term effectively lowers the order 
in to which this term contributes (naively, to order ^^^'^)- This is, however, not the case since the 
leading contribution is canceled by a similar one arising from the first term in Eq. (j20p . Therefore, 
Eq. (jl9p generates all terms in the expansion of the integral Ii up to and including 0{(,'^). 

For the sake of completeness, we list below terms in the expansion of /(p,l) in powers of ^, which 
are responsible for the contributions of orders ^^/^ and to the expansion of Ii: 

h + f^-h = ••• + e^^(P,i)Q-^-^ + ^)+C^<i>.(p)(^ + ^)+e^^4(p)^ 

e^cD,(p)b4 g2,^^(p)fc6 ^3^^(p)fc6 

2(12 + ^52) p + ^52 + 2(12 + ^^2)2 ' i I 

where ellipses refer to terms already given in Eq. (jl9p and the quantity $4(p) is defined as 

^4(p) = ^^V?VfcI>(p,l)|,^„. (22) 

Higher-order contributions in ^ can be systematically obtained along the same lines. In particular, 
one verifies that the low-1 regime yields only terms with half-integer powers of ^, whereas integer 
power terms occur from the high-1 regime. The intermediate region plays a role of regularization 
leading to scale-less integrals which cancel the ultraviolet-divergent (infrared-divergent) terms in the 
low-1 (high-1) regimes. 
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4.2 Results 



We are now in the position to apply the expansion method described in the previous subsection. In 
particular, one immediately observes that there is no recoil correction at order As seen in the 

previous subsection, non-integer powers of ^ can only appear from the expansion of the integrands in 
the low-1 regime, see Eq. (I13p . First, we note that the correction Alst given by Eq. (jlOh does not 
yield the non-integer powers in Thus, performing the expansion in Eq. ([8]), the recoil correction for 
the isovector KN case at order ^^/^ read^ 



A/i = -h = 2cJ-^ (P) - (P)] (P^ + 7^) / 



(2vr)3l2i2 + 2^(p2 + ^2)) (23) 



For the isoscalar case, the integral I^n contributes at the same order as Iq. At order the recoil 
correction A/q can be rewritten as 



A/o=/o+/..=^y^^^(p) 

where G^v^ is the full NN Green function defined as 



G..(p,q,^(l))-§^ 



^ (q) , (24) 



r (r, F\= '^(P - q) , 1 MjVjv(p,q, E) 

L.NN{P,^,^) ^2/^^_E-iO^ A{2Tr)3m%{pymN-E-iO){qymN-E-iO) ^ ' 

i2 

and E(\) = —e^ — 2m]^ this order. On the other hand, the Green function Gnn can be expressed 

in terms of a complete set of the bound- and continuum-state wave functions ^'(p) and ^'^^''(p), 
respectively, which are eigenvectors of the two-nucleon HamiltoniarH 

—Sd — E J k^/mAT — — zO 
The wave functions V'(p) = {^(p); ^i"*" ■*(?)} satisfy the Schrodinger equation 

(■AT 



L _ ^(p) = _ J_ f p^Vip, p')^(p') , (27) 



with e = {— e^^; fc^/mTv}, respectively. The normalization for the potential is chosen so that the 
Lippmann-Schwinger equation for Mat at is given by 

MNN{p,ci,E)-V{p,<i)-^J (2,)3 E-p'ym,, + zO ' ^'^^ 



Substituting Eq. (1261) into Eq. (|24|) one obtains 

, dHd^k 1 



(27r)6 l2 + 2^(k2+72) 



I d3p^(p)xI/(+)t(p) 



0, (29) 



*The integral over 1 yields a contribution ~ leading to the correction for /o,i at order 

^Note that the high-1 regime is perturbative, whereas the low-1 regime is not. Under this we mean that, e.g., in the 
low-1 regime one has to consider the full A'^A'^ amplitude in Eq. (|25p without expanding it in Born series - since all terms 
in this expansion contribute at the same order in On the other hand, if we are extracting integer powers of ^ in the 
high-1 regime, the perturbative treatment of the nucleon-nucleon scattering amplitude is legitimate, see Eq. (|30|l . 
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Figure 2: Recoil corrections in the double-scattering process for the isovector (left panel) and isoscalar 
(right panel) cases. The solid curves correspond to the results of full numerical calculations (see text 
for more details), the dashed curves represent the first non- vanishing recoil corrections at order ^. The 
arrows indicate the results for -rrd- and ^d-scattering. 



where we exploited the othogonality of the bound- and continuum-state wave functions and made use 
of the normalization of ^'(p) according to / (i^p^^(p) = {2tt)^. Thus, at order ^^/^ there is a complete 
cancellation of the recoil corrections both for isoscalar and isovector types of KN interaction. The 
origin of the cancellation for the isovector case is explained by the Pauli principle (see Ref. |26j). 
whereas for the isoscalar case it is the orthogonality of the bound-state (deuteron) and continuum- 
state (A^A^ intermediate state) wave functions in the '^Si partial wave. Note that for vrd-scattering 
similar cancellations were observed in Ref. [25] using a potential-model approach. 

The non-vanishing recoil corrections to the static term appear at order ^ both for isoscalar and 
isovector KN interactions. In order to perform an analytic study of this and higher-order terms 
in the expansion we choose NN interaction in the separable form V{p,p') = ^gip)g{p'), where 
5(P) = (p' + P^)-\ ^(P) = Ng{p){p^ + 72)-\ N = V87r7/3(7 + P)^ and /? = 1.4 fm"!. 

Performing the Fourier transform and making use of the Schrodinger equation, one obtains the 
following result for the linear corrections: 

A4 = — [ d^rr^Mh'^ - A)^(r) [ d^rd^r'^(r)^(r')Vir,r')\r - r'\ . (30) 

4tt J IQirmN J 

Evaluating these terms for the employed A^A^ interaction we get 

A/f A/I 

— J-«0.6^ and —^^-0.3^. (31) 

1st 1st 



In Fig. [2] we show the results for the recoil correction (in units of the static term 1st) for isovector 
(left panel) and isoscalar (right panel) contributions as a function of ^. The results of our full numerical 
calculation of the recoil corrections, see Eqs. ([9]) and (llOp . without expanding in ^ are shown by 
the solid curves. Surprisingly, even for Kd scattering the nucleon recoil effect turns out to be not that 
large as one could a priori guess. As can be seen from the figure, the nucleon recoil for the double- 
scattering process amounts just to 10-15% of the static contribution. To understand the origin of the 



11 




Figure 3: Convergence of the expansion in ^ for the recoil corrections in the double-scattering process 
for the isovector KN interaction. The notation is explained in the text. The arrows indicate the 
results for ird- and Xd-scattering. Note that the (trivial) kinematical pre- factor 1/(1 + C) was not 
expanded in ^ (see the text for details). 



smallness of the effect we compare the full result with the EFT calculation based on the expansion in 
powers of ^. Our results at order ^ are shown by the dashed lines. As visualized in Fig. [21 the full 
recoil correction changes its sign in the considered interval of for the isovector case. Here, the linear 
approximation fails completely to describe the full result. On the other hand, for the isoscalar case 
the recoil correction has a constant sign, and the linear approximation yields a reasonable result. 

To further explore the convergence of the expansion in ^ in the case of the isovector KN interaction 
we calculated higher-order corrections using the expansion method described above. The results 
are shown in Fig. [3l The thick red line corresponds to the full result, solid lines represent the 
calculations up to and including half-integer powers of whereas the results shown by dashed lines 
include, in addition, the next-higher integer power of ^. Fig. [3] demonstrates that already at order 

one reproduces the bulk of the effect while the order-^^ calculation already provides a very good 
approximation to the underlying result for Kd scattering. One can see from this figure that the 
smallness of the net recoil effect is accounted for by the specific cancellation pattern among different 
recoil corrections. In particular, there is a huge cancellation between the first integer (at order ^) and 
the first non- integer (at ^^/^) corrections that even leads to the change of sign for the recoil effect. 
Further, while improving convergence at smaller ^, the inclusion of higher-order half- integer terms 
results in an oscillatory behavior around the full result at larger ^. Actually, looking at Eqs. (jl9p and 
(j21|) . one can already see that the sign in front of the leading non- integer terms changes while going 
from order ^^/^ to ^^^^ (cf. the last term in Eq. (I19p and the last three terms in Eq. (j2ip ). This can be 
explained by considering the expansion of the 3-body propagator in Eq. (jlip in the low-momentum 
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Figure 4: Inclusion of the recoil corrections in the multiple-scattering process. 



regime, which leads to an alternating series 

/ = (l^*. + l% + ..0^^(l-^p^ + j(p^^-...). (32) 

Due to this pattern, terms with half-integer (HI) powers of contributing to /i also show alternating 
behavior 

II™ = -^{2.lf/^ - 0.96^^/2 + 0.85^^/2 _ o_sif/2 + 0.8^"/^ - 0.8^^^/^ + ...), (33) 
1 + ? 

which results in their partial cancellation. Notice further that the (trivial) pre- factor 1/(1 + ^) in Ii, 
see Eqs. (Ilip and ([33]), produces large coefficients when expanded in ^. After expanding it in powers 
of ^, Eq. ([33]) turns into 

/i™ = (2.1^3/2 _ 3.07^5/2 _^ 3.91^7/2 _ 4 72^9/2 _^ 5,53^11/2 _ 6.33^^^/^ + ...). (34) 

The convergence in the latter case is much slower (albeit the series still converges at the value of ^ cor- 
responding to the kaon mass). For this reason, we kept this trivial kinematical pre- factor unexpanded 
when showing our results in Fig. [3l 

Finally, we would like to emphasize that due to large cancellations between individually sizable terms 
discussed in this paragraph, the recoil effect turns out to be not that large even for Kd scattering. 
We found that it is about 10-15% of the static piece for the separable model of A^A^ interaction (see 
arrows that indicate the results for Kd scattering in Fig. [2]). 



5 Recoil effect in multiple scattering 

Due to the strong ATA interaction, the recoil effect in the multiple-scattering diagrams might be 
significant and must be studied quantitatively. Moreover, it can be shown that, starting from the 
quadruple ^A scattering process, the recoil correction becomes nonzero already at order ^/^. A 
detailed discussion of the recoil effect in the multiple scattering will be reported elsewhere [30]. Here 
we just sketch the method. The crucial assumption of the method is that the recoil corrections can be 
treated perturbatively, even if the static KN interactions have to be resummed to all orders. Thus, 
one may study an insertion of 1,2,. . . "retarded" blocks in the diagrams that contain infinitely many 
static kaon propagators. An example is shown in Fig. H] where we have dressed the double-scattering 
block - studied in detail in the previous section - by static kaon rescattering in the initial and final 
states. In order to keep track of various powers of ^ within this method, it is essential to have an 
explicit perturbative expansion of the retarded block that can be achieved by using the technique 
described in the present article. 
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6 Conclusions 



We have studied the nucleon recoil effect for Kd scattering using EFT. Specifically, using the expansion 
method of the Feynman diagrams in EFT, we have calculated recoil corrections to the double-scattering 
process in a systematic expansion in the half-integer powers of the parameter ^ = Mk/itin- It was 
shown that the leading correction to the static term, which emerges at order cancels completely 
both for isoscalar and isovector types of KN interaction. The origin of the cancellation for the 
isovector case can be explained by the Pauli principle whereas for the isoscalar case it stems from the 
orthogonality of the bound state (deuteron) and continuum (A'^A^ intermediate state) wave functions 
in the partial wave. The coefficients of higher order terms in the expansion in ^^/^ appear to be 
of a natural size and the series converges even for the value of corresponding to the physical kaon 
mass. Due to a significant cancellation that takes place between individually sizable terms at different 
orders, the net recoil effect for the double-scattering process is found to be just about 10-15% of the 
static contribution for the separable model of NN interaction. We also sketched the method that can 
be used to include the nucleon recoil for the multiple-scattering process. A more detailed discussion 
of this and other aspects will be presented in the forthcoming publication [30|- 
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